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, Entanglement of atomic systems is a promising resource for performing various computational and communication 
tasks originated in quantum information theory. Recent experimental progress has been achieved in the field of 
CO ' creation of multipartite highly nonclassical states in a linear Paul trap [l|, . Production of such states requires a 
^ . nontrivial implementation of a sequence of nonlocal operations (unitary gates) that entangle originally uncorrelated 
atoms in a desired way. Apart from these sophisticated methods for manipulating the atomic states with the help of 
a precisely tuned laser field, it has been predicted that originally independent atoms may become entangled after the 
■ interaction with a thermal electromagnetic field. At first sight this may seem counterintuitive because thermal light is 
', considered to be a decoherence source. But indeed, if we stay in the quite modest case of entangling just two atoms, 
\0 • Kim et al. Q as well as Bose et al. [1] and An et al. Q have reported the occurrence of entanglement of initially 
\ separable mixed states after the interaction with thermal light. These studies belong to a broader group of works on 
. the behavior of an atomic ensemble coupled to different kinds of bosonic environment 0, 0] . 

Loosely speaking, all previous studies of the coupling of an atomic system to a thermal bath 0, [^, Q concluded 
^ ' , , with the interesting result that entanglement not only does not disappear after the interaction but it can be even 
^ ' created from originally separable states. In this work we go further. In the framework of the model describing 
' the interaction of an atomic ensemble and a single mode photon field worked out by Tavis et al. Q we study two 
^ • physically relevant families of initial pure atomic states coupled to a thermal bath. The states can be tuned by 
changing its Schmidt number (characterizing the degree of entanglement of initial states) from factorized states to all 
four maximally entangled states (Bell states). The calculations give quite unexpected results that, to our knowledge, 
have not been reported so far. It is found that for an input state of the form \^ = sin/3 |00) -I- cos/3 |11) and 
for a sufficiently 'hot' environment, there is a stable nonzero probability of producing a maximally entangled state 
' irrespective of the value of /3 and, what is more interesting, despite of the fact that the resulting state after the 
interaction is separable. This procedure is usually called postselection. Thus we present here a probabilistic source 
of maximally entangled atomic pairs. 

A detailed study of the family of input atomic states |$^) = sin/3 |01) + cos/3 |10) shows that entanglement of output 
states is higher in comparison with input states for a wide interval of (3 values. The picture started in Refs. [3, Q is 
now completed by showing that the entangling properties of a thermal field prevail even at very low temperatures, so 
that vacuum fluctuations can generate partially entangled states from initially factorized states. 

The structure of the paper is as follows. In Section |T] we introduce the entanglement of formation (EOF) [l^. This 
entanglement measure is calculated for general two-qubit states with the help of the analytical expression for the 
concurrence In Section [IT] we describe the Tavis- Jaynes-Cummings model setting the stage for the main results 
presented in two subsections where the evolution of thermal light with two different atomic input configurations |$^) 
and \^ 13) is studied. The details about the calculations are given in Appendix. 
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I. ENTANGLEMENT MEASURES 

Most axioms required for entanglement measures [l^l are accomplished by the entanglement of formation (EOF) 
defined by 

Ef{q)^ inf Y.P^EiM^U), (f) 

where E {•) ~ S{Tti(^2) [|</')(</'|i]) is called the entropy of entanglement (the von Neumann entropy of the state of 
interest that is traced over one of its subsystems). The EOF is generally uneasy to calculate even for the lowest 
dimensional systems. However, Hill and Wootters proved [ll| that the EOF can be calculated analytically in the case 
of an arbitrary two qubit state g. They showed that 

EFig) = h(^^ [l+^l-C^ig))y (2) 

where h(x) — — a;log2(a;) — (1 — x) log2(l — x) is the Shannon entropy and 

C{g) = max{0, Ai - A2 - A3 - A4}. (3) 

The parameters A.; in Eq. ([3]) are eigenvalues in descending order of the square root of the matrix 

R{g) = g{ay®ay) g* {ay ® ay) (4) 

with * representing the complex conjugation in the standard basis and ay is the y-Pauli matrix. C is called the 
concurrence and since the EOF is monotonous in C, the concurrence can also be considered as the entanglement 
measure. Nevertheless, as emphasized in [l3| . the proper entanglement measure is the EOF. Using the concurrence 
for this purpose could cause confusion when investigating entanglement of systems with dimensionality nxn for n ^ 2 
since in that case the definition of concurrence is not unique. 

Also in some cases, e.g. in Q, the negativity [3l as an entanglement measure is preferred. The reason is that it is 
easier to calculate (tracing over a partially transposed density matrix). On the other hand, the negativity does not 
coincide with the entropy of entanglement introduced above for pure bipartite states where it is the unique measure 
of entanglement [l^. Rather recently [TB], connection of the negativity with the concurrence (and thus to the EOF 
that coincides with the entropy of entanglement for pure bipartite states) was discovered in terms of convex roof 
construction. Anyhow, in this paper the particular forms of the density matrices allow us to calculate the EOF 
directly from the concurrence. 

The second note concerns the term Schmidt number. It characterizes the degree of entanglement of input states 
but does not uniquely determine the behavior after the interaction. In other words, even if two input states have the 
same Schmidt number, for example two Bell states, their evolution and thus entanglement after the interaction may 
be completely different. 



II. INTERACTION OF TWO ATOMS WITH THERMAL LIGHT 

The Hamiltonian of two identical two-level atoms interacting with a single-mode electromagnetic field in the dipole 
approximation and standard notation is given by 

H = Hq + Hint 

= nn (b^b + + f (a?^ + + ng -fb + a«5t^ , (5) 

where resonance of the photon energy and the atomic level splitting is assumed, and the rotating-wave approximation 
is used. The analogous Hamiltonian for A^-atoms can be analytically diagonalized as already shown in Ref. 0. Here 
we use the dressed states formalism to obtain the eigenvalues and eigenvectors. Since [Ho,Hint] = 0, Hint induces 
transitions only between the degenerate states of Hq that, for a given field excitation number n, constitutes the tetrad 
{\n) |11), \n + 1) |10), \n + 1) |01), \n + 2) |00)}„. The bare states form a 'semilogical' basis where, e.g., \n + 1) |10) 
represents the state with (n -I- 1) photons, the first atom in the excited state (logical state one) and the second one 
in the ground state (logical state zero). Expressing the interaction Hamiltonian in this basis we get a block diagonal 
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matrix with the n-th block given by 
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After normahzation of the eigenvectors of this matrix, the dressed basis states read 
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In particular, the atomic state l/\/2(|10) — |01)) is a dark state. 
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A. Initial Atomic State |"I>,3) 



In this section, the system light-atoms is assumed to be initially in a factorized state > (/3, 0) = Qught ® Qa 
where the atomic state is a pure state 

1$^) =sin/3|01)+cos/3|10) 
with the Schmidt number sin 13 and the light considered to be in a thermal state 
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with n = [ef'^/kT _ i] the average number of photons. After a direct but otherwise tedious calculation the time 
evolved density matrix p*^*) (/?, t) can be obtained. Since we are interested in the entanglement behavior of the atom- 
atom system, we trace over the electromagnetic field states. The resulting time-dependent reduced density matrix is 
of the form 



^W(/3,t) 



/ai 
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(10) 



The specific expressions of the matrix elements can be found in Appendix. 

The structure of is very friendly for the calculation of the concurrence with Eq. ([3]) and consequently the EOF. 
The result is illustrated in Fig. [T] where oscillations of the EOF are noticeable for all (3. There is a /3o <C 1 that 
depends on n such that for the subinterval f3o < P < tt /2 — Pq the minimum value of the EOF is zero and the local 
maxima (revivals) are always smaller than the initial value for n > and equal for n = 0. Thus, the EOF of the 
atomic entangled states is reduced in the presence of thermal light. As the temperature increases the temporal width 
of the revivals decreases. A brief inspection of the picture reveals other interesting aspects of interacting light with 
this family of states |$^) in the second interval of (3, tt/2 < (3 < tt . Due to the fact that l/-\/2(|10) — |01)) is a 
dark state the behavior is completely different. The EOF exhibits very soft oscillations and is never smaller than its 
initial value. We thus see that in this region, the resulting state is stable under the interaction with thermal light. 
This could be an important hint for quantum engineers (and protocol designers generally) telling them that there 
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(a) (b) 

FIG. 1: The EOF of the output state given by Eq. (|10|l for interacting atomic states |$^) with; (a) a zero photon field, and 
(b) a thermal field with a nonzero average photon number (n ~ 0.64) . A completely different behavior for < /3 < 7r/2 and 
''"/S < P < IT can be seen. EOF revivals are observed in the first interval where, generally, entanglement of input states is 
devalued except for initially factorized or wealily entangled states. On the other hand, in the second interval the EOF is never 
lower than its initial value. We can observe significant values of the EOF for initially factorized atomic states and n = 0. 



exists an important class of initial atomic states where the interaction with a thermal field means not only no loss of 
entanglement but even its enhancement without strong fluctuations. The works presented so far on this topic (e.g. 
Refs. [1, 0]) were limited only to input atomic states with small values of the Schmidt number < /3 < /3o ^ 1 and 
recognized the presence of the dark state but do not discuss the stability of entanglement in the whole range of (3. 
For n = the matrix elements of Eq. pop have a particulary simple form 
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By examining the density matrix for /3 = we may exactly see the influence of the vacuum fluctuations on a factorized 
input state. We see that there is a 50% probability of finding the system with both atoms in their ground state and 
a single photon, while there is a 50% probability that this photon has been reabsorbed by one of the atoms leaving 
the atomic system in the entangled state. As a consequence, spontaneous decay may be regarded as a causative 
mechanism for inducing entanglement between the two two-level atoms when their initial state |$^) is practically 
unentangled (/3 ~ 0) or to enhance it for different /3 in the interval 7r/2 < /3 < tt. 
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FIG. 2: The continuous lines in both plots are the EOF for mixtures of the output states occurring within the time interval 
where the EOF is periodic (n = 0). The dash-dotted lines correspond to the EOF for t = (initial states). In the upper plot 
the interval /3o < /3 < 7i"/2 — /3o with /3o = 1/2 arcsin(l/7) demarcates the region with lower entanglement in comparison to 
corresponding input states. In the bottom plot we can see that entanglement is not present in the interval (— /9i,/3i), with 
/3i = arctan(l/8), hence disappearing even for initially weakly entangled states. 



The fact that EOF exhibits strong oscillations as already illustrated may limit even its evaluation. Evaluation of 
the entanglement requires knowledge of the matrix elements and these cannot be determined just from one copy of 
an unknown quantum state. One possibility is to have many copies of the same unknown state and perform quantum 
state tomography. Recent findings show that the tomography is not necessary [13, [13 but still does hold that it is 
not possible to measure the entanglement from just one copy of a state. Besides, if the average number of photons is 
so large that the time window for picking the desired state is too long compared to the length of the oscillations, the 
measurement could be interpreted as a random pick from the set of all possible outputs and a proper time average 
must be performed. 

Let us study this question for the electromagnetic field initially in the vacuum state as an illustration. In this simple 
case, the natural period is determined by the commeasurable Rabi frequencies \/%g and \/2g 



/2(l + sin2/3) 0\ 

3-sin2/3 3sin2/3- 1 

3sin2/3-l 3-sin2/? 

V Oy 



(12) 



The resulting EOF can be seen in the upper plot of Fig. [2] and, we are able to confirm the qualitative conclusions 
about the entanglement behavior based on Fig. [T]for the averaged EOF. As mentioned above, we can observe nonzero 
entanglement for initially factorized states. For entangled states with /3o < /? < 7r/2 — /3o the averaged EOF is lower 
than for the input states. In general, for the remaining (3 the EOF average is higher than the initial value of the EOF. 

The extension of this averaging process for the case of an interaction with thermal light with n > requires the 
identificatio n of a na tural time scale since in this case an infinite number of the incommensurable Rabi frequencies 
an-i = ^g\/n +1/2 determines the evolution of the system. One such scale could be found by the following procedure. 
Writing the summations in the expression of g^*'' (/?, €) in an integral form using the Abel-Plana formula, it can be 
shown that the adimensional parameter gt^j2kT /Ml is especially important to understand the evolution of the system. 
In fact, it can be used to define a natural time tq ~ g^^ ^JhVl/2kT to distinguish short and long time effects of the 
interaction (see Appendix for more details). If the average procedure from above is applied, tq could also be used for 
this purpose. 
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1.57 1.57 

(a) (b) 

FIG. 3: The EOF of output state (|14|l for interacting atomic states l^'a) with (a) a zero photon field, (b) a thermal field with 
n ~ 0.64. Notice the small enhancement of the EOF for /3 ~ ■k/2 where n was chosen to provide a maximum value of the 
enhancement. 



B. Initial Atomic State |*^) 



Following the same reasoning as in the previous subsection, but for an initial state of the form 

= sin/3|00) +cos/3|ll) 

interacting with light in a thermal state ([9]), we arrive to a reduced density matrix with the structure 



fhi &2\ 

63 &3 

63 &3 

\b* hJ 



(13) 



(14) 



The particular expressions for the matrix elements are given in Appendix. As illustrated in Fig. [3l in this case a 
resonant thermal electromagnetic field is not able to induce, in general, atomic entanglement for input factorized 
atomic states, nor to enhance it for initial already entangled states. An exception is provided for the states with 
f} ~ 7r/2. There, a very slight EOF enhancement is observed for < n < 2 . This effect is discussed at length in 
Ref. [31 and present, but barely observable, in Fig. [3l For the other values of (3, at most, partially entangled input 
states exhibit revivals of the EOF. Thus, in comparison to the input state |$/3), the entanglement properties of a 
thermal field are not so striking. This qualitative observation is confirmed for n = by the mixture calculation in the 
spirit of the previous section. It gives us 



12 



/4(2-cos2/3) 



V 4 sin 2/3 





cos^ /? 
cos^ 13 



4sin2/3N 
cos^ (3 
cos2 (3 

6 cos^ / 



(15) 



7 




P [rad] 



FIG. 4: Depicted values of parameter pi for n « 5.9 which determines the probability of the successful postselection of a 
maximally entangled state l/\/2(|10) + 101)) from outgoing state (|14[l . It is possible to postselect the state irrespective of /3 
and thus it can be done even if the output state is separable. 



The resulting EOF is depicted in Fig. [5] (bottom plot). We see that in the mixture sense, there is not only a separable 
output for initially unentangled states, as discussed above, but also for some of the input partially entangled states. 
Also, if we compare the EOF for the initial state (dash dotted line) and evolved state (continuous line) the output 
states are much less entangled in general. 

Nevertheless, returning to the evolved state another interesting aspect appears. Eq. can be written as 

^(*)(/3,^)=Pl^?l + (l-Pl)^?2 (16) 

where pi = 63 and 

gi = i|01 + 10)(01 + 10| 

92 = -47-(^i|00)(00| + 64|ll)(ll| + 62|00)(lli + 6*|ll)(00|) (17) 

Oi + 64 

and both states can be unambiguously discriminated. To show this, let us consider a projectively measuring apparatus 
with three possible outputs. Two of them project into orthogonal subspaces spanned by |00)(00| and |11)(11| and the 
third one is the complement. Then, from the form of p6|) it follows that if we do not detect anything on the first 
and the second output the state is in the required maximally entangled state |01 + 10)- The probability of this 
procedure is given by the parameter pi. Thus, no additional quantum resource for the discrimination is needed [l9| . 

Therefore, even if the EOF is zero (i.e. when an output state is separable) the entangled pair can be extracted in 
terms of postselection. Also, it is interesting that the maximum value of pi is almost the same not only for different 
input states (see Fig. 2]) but also for different average photon numbers of thermal light. It is worth remarking that 
the higher photon average number is used the less fluctuating the parameter pi we get. This constant character of pi 
is broken only for the case of a zero photon field where, not surprisingly, for an input state |00) we cannot extract a 
maximally entangled state by any means, see Eq. p5|) . 

Note that a similar procedure cannot be applied to the resulting state in Eq. pop from the previous subsection. 
There, we are not able to find a decomposition of a density matrix in the form g ~ gi (B g2 SuS in the current case, 
which ensures that if two arbitrary density matrices are positive the third must be so. 
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III. CONCLUSIONS 



In this work, we have explored the possibihty of preserving or enhancing entanglement between two atoms after 
their interaction with a thermal field. We have investigated the interaction for two kinds of initial atomic states 
|$/3) = l/V2(sin/?|01) + cos/3|10)) and I*;?) = l/v/2(sin/3 |00) + cos/3 Although this physical system had 

been discussed previously, new interesting phenomena have been found and the already reported results have been 
confirmed. For the initial state |$/3), wc observe a very different behavior after the interaction with the thermal 
field for two distinguished intervals of input states. Qualitatively, for almost the whole range (3 = (0,7r/2), the EOF 
is decreased while for the remaining input states there are oscillations in the EOF so that they are never smaller 
than their initial value. This result is more significant at low temperatures, hence, the entanglement enhancement 
mechanism is induced by spontaneous decay. 

For the second family of input states parameterized by \'^fj), there is no entanglement enhancement for almost any 
/3. In this case, the entanglement is generally degraded matching our intuitive conception of an interaction between a 
correlated pairs of qubits and thermal light. Nevertheless, the output atomic density matrix can be decomposed into 
10) = l/\/2 |01 + 10) and a partially entangled mixed state, both of them living in mutually orthogonal subspaces. 
Thus they can be discriminated with zero probability of error. If we consider that this holds for all /3 when fi 
is sufficiently high, it follows that we are able to extract the maximally entangled state |0) even if the input is a 
factorized state and the output is in a separable form. We may conclude that we have at disposal a probabilistic 
source of maximally entangled atomic pairs, when our deposit is just a sufficiently hot thermal light with two atoms 
in a factorized state. 

In current experiments with trapped ions, the energies involved in the transitions compared with room temperatures 
are such that Ml/kT <C 1 for the Innsbruck group (729 nm atomic transition in '*°Ca+) Q and hfl/kT ^ 1 for the 
NIST group {^Be+ at 1.2 GHz) [l|. Thus, one would expect that the former is a good scenario for testing the 
predictions made for the initial state l^^g), in particular, the entanglement enhancement. Meanwhile, the NIST group 
parameters could be appropriate to study the possibility of stabilized postselection of the maximally entangled state 
described for the interaction of thermal resonant light with l^'/j). 
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APPENDIX 
Initial Atomic State |<1'^) 

With the help of the inverse transformation from the bare state to the dressed state basis, we rewrite the initial 
density matrix determined by the atomic state |$^) and a thermal single-mode field 

oo 

^.W(/3,t = 0) = Et[(^ + '"'^'^)("'^>" + '^>")("^^'" + ^^'") 

~a/2cos2/3 ((- |3)„ + |4) J (2 |„ + |2>„(~ (3 |„ + (4 |J) - 2(1 ~ sin2/3)|2)(2|„] , (A.l) 

where p„ = [1 — exp {—hn/kT)] exp {~nMl/kT). The time evolution of the whole system is directly found. Tracing 
over the resulting density matrix we obtain an output atomic state (Eqs. (jlOp l with elements 
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with 



= — ^^^^cosa„_ii 
— i n + 1/ 2 



n=0 
oo 



h2{t) = y^p„ COSCtn-it = -^h'l(t). 



(A.7a) 
(A.7b) 



ttn-i = 2gy^n + 1/2 is the Rabi frequency coming from eigenvalues of the interaction Hamiltonian ([6]) and rn± 
1 ± Marctan (exp {-Ml/2kT)) where M = [1 - exp (-?ifi/fcT)] exp {hn/2kT). 

Expressions (jA.7[) appearing in the matrix elements can be written in a more compact form 

oo oo oo 

h2{t) = 2 sinh k ^ e''(^"+^^ cos V2n + Igt = 2 sinh 2k ^ e"" cos ^gt - ^ e^*^" cos V2ngt (A.. 



n=0 
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n=0 



where k = Using the Abel-Plana formula 



f: F{n) = r F(x)dx + \f{0) + ^ H EMzItl^la, 



n=0 



(A.9) 



the following integral expression for h2{t) is found 

sinh K 
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with an adimensional variable i = that establishes a natural time scale tq = ^ for the description of the 
interaction. Notice that the second integral is an exponential decreasing function of the temperature so that for a 
sufficiently hot environment 



, , , „sinhK _^ - ^sinliK 
h2[t) ^ 2 / e cosiVa^da; = 2 
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The latter integral can be recognized as an error function with an imaginary argument. An integral expression for 
hi{t) can be found in a similar manner. 



Initial Atomic State j^JZ/j) 
Calculations similar to the case of the initial state |$^) lead to the following matrix coefficients of Eq. 
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